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Quantum Anomalous Hall Effect in Kagome Ice 
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An anomalous Hall insulator without magnetic long-range ordering is theoretically reported in 
the absence of the relativistic spin-orbit coupling. It is realized in itinerant electrons coupled with 
the Ising spins on a (111) kagome plane of pyrochlore spin ice in applied magnetic field. We find 
that the kagome-ice type local spin correlation in the magnetization plateau state opens a charge 
gap without magnetic ordering, which results in quantization of the Hall conductivity. By Monte 
Carlo simulation, we identify the anomalous Hall insulating region in the magnetic phase diagram, 
in addition to another anomalous Hall insulator in a fully-saturated state. 

PACS numbers: 73.43-f, 71.30.+h, 71.10.Fd 



Competing interactions in geometrically frustrated 
magnets often give rise to peculiar local spin textures 
or objects composed of several spins. A representative 
example is the two-in two-out spin configuration in each 
tetrahedron in spin-ice compounds [l|, Q . The local tex- 
tures or objects have been a matter of intense studies 
as they are the source of interesting properties in frus- 
trated magnets, such as macroscopic degeneracy with 
residual entropy and a characteristic power-law spin cor- 
relation [1,13. 

Recently, metallic pyrochlore oxides have opened up a 
new aspect in the study of geometrically frustrated mag- 
nets. In these compounds, localized spin-ice moments 
on a pyrochlore lattice interact with itinerant electrons; 
hence, it is expected that the local spin objects affect 
the electronic and transport properties. Indeed, various 
peculiar transport properties were reported, such as un- 
conventional anomalous Hall effects in Nd2Mo207 [5[ and 
Pr2lr207 [sl-lsl, and resistivity minimum in Pr2lr207 [ol 
and Nd2lr207 [lo|. Effects of local spin objects on these 
phenomena were also studied theoretically [ll|, [l2[ . 

The interest in the interplay between local objects and 
itinerant electrons has spanned to a wider range of sys- 
tems. For instance, extended Falicov-Kimball models on 
frustrated lattices show non-Fermi liquid behavior [l3| 
and peculiar metal-insulator transition [ij]. Another re- 
lated study was done on a frustrated double-exchange 
(DE) model, focusing on the nature of loops emergent 
from local spin textures [l5| . These results indicate that 
such interplay offers a fertile ground for exploring uncon- 
ventional electronic and transport properties. 

In this Letter, we propose that a peculiar spin texture 
gives rise to an anomalous Hall insulator (AHI) without 
a magnetic long-range order. We show that the so-called 
kagome ice, which appears in a (111) kagome plane of 
the pyrochlore spin ice in applied magnetic field, affects 
the electronic state significantly by opening a charge gap. 
We numerically show that the insulating kagome ice state 
exhibits a nonzero quantization of the Hall conductivity 
induced by the spin scalar chirality, indicating that it is 
an AHI characterized by the first Chern number. Increas- 



FIG. 1. (color online), (a) Schematic picture of spin ice on a 
pyrochlore lattice. The shaded plane indicates a (111) kagome 
plane, (b) One-in two-out and (c) all-out type spin configu- 
ration. 



ing the magnetic field, we find that the AHI state turns 
into another AHI with saturation of the magnetization; 
this can be viewed as a topological change between the 
two insulators. The conventional quantum Hall state is 
associated with the formation of quantized Landau levels 
in external magnetic field. Recent studies on topological 
states of matter revealed that a quantum (spin) anoma- 
lous Hall effect occurs in the presence of the relativistic 
spin-orbit coupling 16-20|. Besides, a noncoplanar mag- 



netic order can also give rise to an AHI [2l|, |22|. Our 
present result offers yet another AHI without Landau 
levels, spin-orbit coupling, nor magnetic ordering. 

We consider spin-ice type Ising moments on a kagome 
lattice. The localized Ising moments correspond to those 
on a (111) kagome plane of the spin ice model, as shown 
in Fig.lIJa); the anisotropy axis of the Ising spin is along 
the direction connecting the centers of neighboring two 
tetrahedra. When the Ising spins interact with each other 
by the nearest-neighbor (NN) ferromagnetic interaction, 
the ground state remains disordered with macroscopic de- 
generacy, in which all the triangles obey one-in two-out 
or two-in one-out spin configurations [Fig. [TJb)]. An ex- 
ternal magnetic field perpendicular to the plane partially 
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lifts the degeneracy by enforcing the upward (downward) 
triangles to be two-in one-out (one-in two-out), but the 
ground state is still disordered. This is called the kagome 
ice, which is realized in the magnetization plateau state 
in spin ice (23l-[25|. For a higher field, all-in(out) up- 
ward (downward) triangles (magnetic monopoles in the 
pyrochlore spin ice [26|) are introduced into the kagome 
ice, resulting in a first-order transition to a saturated 
state. 

To clarify the effects of such characteristic local spin 
textures on itinerant electrons, we here consider a single- 
band DE model 13, [HI with the Ising spins on a kagome 
lattice; 

H = - ^{Ujc\cj +H.C.) - ^hzCosOi. (1) 

Here, q {c\) is the annihilation (creation) operator of 
an itinerant electron at iih site, whose spin index is 
dropped as the spin is completely aligned parallel to 
the localized spin at each site. The anisotropy 
axis of the localized spin depends on the sublattice; 
Si = {Sf ^ Sf , S^) = S {sin 6i COS (fi^ sin 6i sin (pi^ COS 6i), 
where (Oi^ipi) = (arccos(±|), ± f) for the sub- 
lattice Us = 1,2,3 and 6' = 1 [see Fig. [U^c)]. The ef- 
fective transfer integral Uj depends on the relative an- 
gle of neighboring Ising spins, which is given by tij = 



t(cos - 



■ cos - 



^). The sum is 



taken over NN sites on the kagome lattice. In Eq. ([T]), 
the effect of external magnetic field is taken into account 
by the Zeeman term only for Ising spins for simplicity. 
Hereafter, we set the energy unit t = 1, the length of 
Bravais lattice vector a = 1, the Boltzmann constant 
/cb = 1, and the unit of conductance e'^/h = 1 {e is the 
elementary charge and h is the Planck constant). 

In this model, the spin-charge coupling induces an ef- 
fective ferromagnetic interaction between the localized 
spins at general filling 27, HH. Although the NN fer- 
romagnetic interaction favors macroscopic degeneracy as 
described above, the spin-charge coupling also gives rise 
to further- neighbor interactions, which generally lift the 
degeneracy and select a magnetic order in the ground 
state. However, as the further- neighbor interactions are 
usually much weaker than the NN one, the liquid-like 
state with strong local correlations of two-in one-out or 
one-in two-out types is expected to emerge in the inter- 
mediate temperature (T) region between a high-T disor- 
dered state and the low-T ordered phase. This is indeed 
confirmed by the finite-T Monte Carlo (MC) simulation 
as we discuss later. In the following, we focus on how 
the electronic and transport properties in the liquid-like 
state evolve in applied magnetic field. 

First, for simplicity, we study the problem by taking 
a simple average (arithmetic mean) over different config- 
urations of spins. Instead of the magnetic field /i^, we 
control the net magnetization along the z direction per 



triangle, = where is the number of 

sites. The average is taken over different spin configu- 
rations with a fixed ratio of different types of triangles 
as follows. For < < 1/3, we consider the manifold 
in which all the triangles are of two-in one-out or one-in 
two-out type; their ratio is controlled so as to realize the 
given value of m^. At = 1/3, all the upward (down- 
ward) triangles are in the two-in one-out (one-in two-out) 
configuration, which is the kagome ice. For > 1/3, 
we introduce all-in upward and all-out downward trian- 
gles. Eventually, at the saturation to = 1, the sys- 
tem shows a long-range order of alternating all-in and 
all-out triangles. At each value of < 1, we calculate 
the electronic state of the model in Eq. ([1]) by numeri- 
cal diagonalization with taking the simple average over 
different spi n configurations generated by a loop update 
method [1^. The conductivity is also numerically ob- 
tained by using the Kubo formula (29| . 

Figures [2fa) and Efb) show the electronic density of 
states (DOS) for different values of m^. The solid and 
open arrows on the spectra indicate the Fermi levels for 
the electron filling n = 1/3 and n = 2/3, respectively 
= Figure a) shows the results for 

< ruz < 1/3, where all the triangles are in two-in one- 
out or one-in two-out configurations. The result indicates 
that, although there is no long-range magnetic order, the 
characteristic spatial correlations in the kagome ice mani- 
fold develops a charge gap at n = 2/3 near = 1/3 (soj . 
Figure [2fb) shows the results for 1/3 < < 1, where 
the all-in upward or all-out downward triangles (mag- 
netic monopoles) are introduced into the kagome ice man- 
ifold. At the Fermi level for n = 1/3, the introduction 
of monopoles leads to a dip in the spectra, and opens a 
full gap for rriz ^ 0.7. On the other hand, the energy 
gap at n = 2/3 and = 1/3 is closed by introducing 
monopoles. However, for rriz ^ 0.7, DOS develops a dip 
again and shows an energy gap at 1 [sij . 

To gain further insight into the peculiar changes, we 
calculated the longitudinal (axx) and transverse {(Jxy) 
conductivities. Figure [2fc) shows the results at n = 2/3. 
For < rriz < 1/3, ctxx shows monotonic decrease in 
accordance with the growth of energy gap in DOS. On 
the other hand, axy shows monotonic increase while in- 
creasing rriz , and surprisingly, the numerical results show 
a quantized value at (Jxy = 1 in the gapped state near 
rriz = 1/3 [si^. For rriz > 1/3, (Jxx increases and shows 
a hump at rriz ^ 2/3; finally it decreases to zero in the 
all-in/all-out insulator at = 1. Correspondingly, axy 
decreases from 1 with showing a sign change at rriz 2/3, 
and converges to another quantized value = — 1 at 
rriz = 1- 

The non-monotonic change of the Hall conductivity 



axy is explained by the Berry phase mechanism [2l|, |33|, 



34]. In the previous study by assuming a q = mag- 
netic order [215, shown that the Hall conductance 
depends on the scalar chirality of localized spins in the 
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FIG. 2. (color online), (a), (b) DOS for different vrtz calcu- 
lated by the simple average. See the text for details. Solid 
and open arrows indicate the Fermi levels for n = 1/3 and 
n — 2/3, respectively. The data for different vrtz are plotted 
with the offset of 0.14 and 0.19 in (a) and (b), respectively, 
(c) shows the longitudinal and transverse conductivities at 
n = 2/3 as a function of m^. 



three-site unit cell, • Sj x S^. A similar mechanism 
was shown to work through the fluctuations [11]. Our 
result of Gxy is understood by considering that two-in 
one-out (one-in two-out) upward (downward) triangles 
bring a negative chirality — each and that magnetic 

monopoles bring a positive chirality +-^^ each. 

The remarkable point in our results is the quantiza- 
tion of Gxy at +1 in the gapped kagome ice state and its 
switching to —1 accompanied by the closing and reopen- 
ing of the energy gap at n = 2/3 and associated hump in 
Gxx' The change in a^y as well as the hump in a^x sug- 
gests a transition between the kagome ice at = 1/3 
and all- in/all-out ordered state at = 1. The latter 
ordered state was shown to be an AHI characterized by 
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FIG. 3. (color online). MC results for (a) vrtz and Xz^ (b) p 
and xp (c) ^, and (d) g^x and a^y at n = 2/3 and T = 0.03. 
The data were calculated for the system sizes ranging from 



the first Chern number C = —1 [21]. On the other hand, 
the kagome ice at = 1/3 is a paramagnetic state with 
constrained local spin configurations; hence, the quanti- 
zation of Gxy as weh as the gap opening in DOS is highly 
nontrivial, suggesting that the liquid-like kagome ice is 
also another AHI. The closing and reopening of the en- 
ergy gap as well as the hump of a^x is interpreted as 
a phase transition between topologically different insula- 
tors with a transient metallic state in between. 

The above calculations suggest that a similar change 
between the kagome ice and all-in/all-out ordered state 
can take place in the model ([1]) in applied magnetic field 
hz. To confirm this, we next investigate the thermody- 
namic behavior by a MC simulation. In this method, the 
MC updates of localized spins are done by directly eval- 
uating the effective action from itinerant electrons using 
exact diagonalization 35, 36| . 

FigureOshows the MC results at T = 0.03 and n = 2/3 
in applied field hz- All the results consistently indicate 



two sharp crossovers while increasing hz ; one is from the 
zero-field state to the kagome ice state by switching on 
hz^ and the other to the all-in/all-out ordered state at 
hz ~ 0.34. Figure [3l^a) depicts the results for and its 
susceptibility Xz- They clearly indicate the presence of 
1/3 magnetization plateau for 0.1 < hz ^ 0.3 followed 
by rapid increase of for larger hz and a peak of Xz at 
hz = 0.341(3). As shown in Fig.[3fb), the ratio of two-in 
one-out and one-in two-out configurations in the system, 
p, stays close to 1 in the plateau regime, while it decreases 
for larger hz] the corresponding susceptibility Xp shows 
a sharp peak at hz = 0.340(3). Figure [3](c) shows the 
net scalar chirality = ^ j /c)(^^ ' ^ ^^)' where 
the sum is taken over all the triangles and the indices 
(i, j, k) are in the counterclockwise order in each triangle. 
It approaches and stays near the minimum value —-^^i 
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pling stabilize some magnetic order, whose ordering pat- 
tern may become complicated depending on hz as well 
as the electron density. It is left for future study to in- 
vestigate the full phase diagram including the low-T or- 
dered phases. When the model is extended to more realis- 
tic situation on the three-dimensional pyrochlore lattice, 
the AHI behavior might be smeared out, but the non- 
monotonic change of axy is likely to remain in applied 
magnetic field. It is interesting to note that such non- 
monotonic behavior of the Hall resistivity was observed 
in the (111) magnetic field in Pr2lr207 [6, 8]. 
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FIG. 4. (color online). Contour plots of the MC results for 
(a) axy and (b) axx in the hz-T plane. The calculations were 
done for N — 3 x 6^ size systems. The white crosses show the 
parameters at which the MC acceptance rate becomes lower 
than 1%. 



followed by rapid increase for a higher field with showing 
a sign change at hz = 0.325(5). Correspondingly, (Jxy 
becomes 1.0 associated with strong suppression of axx in 
the kagome ice plateau region, and shows a sign change 
at hz = 0.328(7) accompanied by a hump of axx- The 
result clearly indicates that the second crossover occurs 
between two topologically different insulating states, the 
kagome ice and all-in/all-out ordered states. 

The two insulating phases are clearly observed in the 
contour plots of axy and axx in the hz-T plane, as shown 
in Fig. Ill Figure m^a) clearly indicates two regimes with a 
finite Hall conductivity axy of an opposite sign: a dome- 
like region centered at hz ~ 0.15, and another one for 
hz ^ 0.34. The two regimes are also clearly indicated 
by strong suppression of axx^ as shown in Fig. ID^b). 
Thus, the phase diagram shows that the model in Eq. ([T]) 
exhibits two different insulating regimes with different 
topological characters, as expected in the analysis by us- 
ing the simple average in Fig. O 

Our comprehensive analyses provide convincing ev- 
idence of the existence of kagome ice AHI in the 
intermediate-T liquid-like regime. This is a peculiar 
topological state of matter in the absence of Landau lev- 
els, spin-orbit coupling, and magnetic ordering. Upon 
further decreasing T, it is expected that effective further- 
neighbor interactions induced by the spin-charge cou- 



[1] M. J. Harris, S. T. Bramwell, D. F. McMorrow, T. Zeiske, 

and K. W. Godfrey, Phys. Rev. Lett. 79, 2554 (1997). 
[2] A. P. Ramirez, A. Hayashi, R. J. Cava, R. Siddharthan, 

and B. S. Shastry, Nature (London) 399, 333 (1999). 
[3] S. V. Isakov, K. Gregor, R. Moessner, and S. L. Sondhi, 

Phys. Rev. Lett. 93, 167204 (2004). 
[4] C. L. Henley, Phys. Rev. B 71, 014424 (2005). 
[5] Y. Taguchi, Y. Oohara, H. Yoshizawa, N. Nagaosa, and 

Y. Tokura, Science 291, 2573 (2001). 
[6] Y. Machida, S. Nakatsuji, Y. Maeno, T. Tayama, T. 

Sakakibara, and S. Onoda, Phys. Rev. Lett. 98, 057203 

(2007). 

[7] Y. Machida, S. Nakatsuji, S. Onoda, T. Tayama, and T. 

Sakakibara, Nature (London) 463, 210 (2009). 
[8] L. Balicas, S. Nakatsuji, Y. Machida, and S. Onoda, 

Phys. Rev. Lett. 106, 217204 (2011). 
[9] S. Nakatsuji, Y. Machida, Y. Maeno, T. Tayama, T. 

Sakakibara, J. van Duijn, L. Balicas, J. N. Millican, R. T. 

Macaluso, and J. Y. Chan, Phys. Rev. Lett. 96, 087204 

(2006). 

[10] M. Sakata, T. Kageyama, K. Shimizu, K. Matsuhira, S. 

Takagi, M. Wakeshima, and Y. Hinatsu, Phys. Rev. B 

83, 041102(R) (2011). 
[11] S. Onoda and N. Nagaosa, Phys. Rev. Lett. 90, 196602 

(2003). 

[12] M. Udagawa, H. Ishizuka, and Y. Motome, Phys. Rev. 
Lett. 108, 066406 (2012). 

[13] M. Udagawa, H. Ishizuka, and Y. Motome, Phys. Rev. 
Lett. 104, 226405 (2010). 

[14] H. Ishizuka, M. Udagawa, and Y. Motome, Phys. Rev. B 
83, 125101 (2011). 

[15] L. D. C. Jaubert, S. Piatecki, M. Haque, and R. Moess- 
ner, Phys. Rev. B 85, 054425 (2012). 

[16] F. D. M. Haldane, Phys. Rev. Lett. 61, 2015 (1988). 

[17] C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 146802 
(2005). 

[18] X. L. Qi, Y. S. Wu, and S. C. Zhang, Phys. Rev. B 74, 



5 



085308 (2006). 

[19] J. E. Moore and L. Balents, Phys. Rev. B 75, 121306(R) 

(2007) . 

[20] L. Fu, C. L. Kane, and E. J. Mele, Phys. Rev. Lett. 98, 
106803 (2007). 

[21] K. Ohgushi, S. Murakami, and N. Nagaosa, Phys. Rev. 

B 62, 6065(R) (2000). 
[22] I. Martin and C. D. Batista, Phys. Rev. Lett. 101, 156402 

(2008) . 

[23] K. Matsuhira, Z. Hiroi, T. Tayama, S. Takagi, and T. 

Sakakibara, J. Phys.: Cond. Matter 14, L559 (2002). 
[24] R. Higashinaka, H. Fukazawa, and Y. Maeno, Phys. Rev. 

B 68, 014415 (2003). 
[25] Z. Hiroi, K. Matsuhira, S. Takagi, T. Tayama, and T. 

Sakakibara, J. Phys. Soc. Jpn. 72, 411 (2003). 
[26] C. Castelnovo, R. Moessner, and S. L. Sondhi, Nature 

451, 42 (2008). 
[27] C. Zener, Phys. Rev. 82, 403 (1951). 
[28] P. W. Anderson and H. Hasegawa, Phys. Rev. 100, 675 

(1955). 

[29] The calculations were done on 4^ superlattices of A/" = 
3 X 18^ sites and T — 10~^ by averaging over 32 realiza- 



tions of different spin configurations. Conductivities are 
calculated with assuming a small inelastic scattering rate 
r"^ = (3 X 18^ X 42)-i. 

[30] The flat band at £ = 1 originates from the bound states 
within the hexagons with all spins "in" or all spins "out" , 
as a consequence of the quantum phase interference. 

[31] In the all- in/ all- out ordered state at rriz — 1, DOS has a 
particle-hole symmetric form with two gaps at n = 1/3 
and n = 2/3 above and below the central flat band [211 ]. 

[32] Gxy becomes +1 for all kagome ice configurations that 
we generated numerically within errors of the order of 
10~^. Small deviations presumably come from a small T 
introduced in the calculations. 

[33] D. Loss and P. M. Goldbart, Phys. Rev. B 45, 13544 
(1992). 

[34] J. Ye, Y. B. Kim, A. J. Millis, B. L Shraiman, P. Majum- 
dar, and Z. Tesanovic, Phys. Rev. Lett. 83, 3737 (1999). 

[35] S. Yunoki, J. Hu, A. L. Malvezzi, A. Moreo, N. Furukawa, 
and E. Dagotto, Phys. Rev. Lett. 80, 845 (1998). 

[36] MC calculations were conducted up to = 3 x 10^ with 
typically 20000 MC measurements after 6000 steps of 
thermalization. Conductivities were calculated by taking 
r = 0.05. 



